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ABSTRACT 

Context. HH 30 is a well-known Pre-Main-Sequence star in Taurus. HST observations have revealed a flared, edge-on disk driving a 

highly-collimated optical jet, making this object a case study for the disk-jet-outflow paradigm. 

Aims. We searched for a molecular outflow, and attempted to better constrain the star and disk parameters. 

Methods. We obtained high angular resolution (~ 1") observations of the dust continuum at 2.7 and 1.3 mm, and of the 12 CO J=2— 1, 
13 CO J=2-l & J=l-0, C I8 J=l-0 emissions around HH30. A standard disk model is used to fit the 13 CO J=2-l ;(v-plane visibil- 
ities and derive the disk properties, and the stellar mass. An ad hoc outflow model is used to reproduce the main properties of the 
12 CO J=2-l emission. 

Results. The rotation vector of the disk points toward the North-Eastern jet. The disk rotation is Keplerian: Using a distance of 140 pc, 
we deduce a mass of 0.45 M Q for the central star. The disk outer radius is 420 AU. A highly asymmetric outflow originates from the 
inner parts of the disk. Only its North-Eastern lobe was detected: it presents to first order a conical morphology with a 30° half opening 
angle and a constant (12 km s~') radial velocity field. Outflow rotation was searched for but not found. The upper limit of the outflow 
rotation velocity is 1 km s -1 at 200 AU of the jet axis. 

Conclusions. HH 30 is a low mass TTauri of spectral type around Ml and age 1 to 4 Myrs, surrounded by a medium size Keplerian 
disk, of mass around 4 10~ 3 M G . It beautifully illustrates the jet-disk-outflow interaction, being so far the only star to display a jet 
and outflow connected to a well defined Keplerian disk, but reveals a surprisingly asymmetric (one-sided) outflow despite a relatively 
symmetric jet. Furthermore, these observations do not enable to assign the origin of the molecular outflow to entrainment by the 
optical jet or to a disk wind. In the latter hypothesis, the lack of rotation would imply an origin in the inner 15 AU of the disk. 

Key words. Individual: HH30 — stars: formation — stars: circumstellar matter — ISM: dust — ISM: molecules radio lines: molec- 
ular 



1. Introduction 

During the last decade, the properties of the circumstellar en- 
vironment of pre-main-sequence (PMS) stars have been exten- 
sively studied. Millimeter arrays, such as the IRAM Plateau 
de Bure interferometer, routinely provide CO line and contin- 
uum images of low-mass PMS stars which show that many of 
them are surrounded b y large (Rant ~ 200 - 800 AU) disks in 
Keplerian rotation (ej > iKoerner et alJll993t Dutrev etaill 19941 
[ Mannings et"afl 119971 iGuilloteau & Dutrevlll998l ISimon et al.l 
2000). For several sources, the angular resolution and sensitiv- 
ity of the observations enable to fit disk models to the data, 
hence constraining physical parameters of the d isk, and even 
the mass of the central source ( iDutrev et all2006l and references 
therein). While mm continuum observations are sensitive to the 
dust thermal emission of the disk, near-infrared (NIR) or opti- 
cal observations can be used to tra ce the stellar light sc attered 
by particles at the disk surface (e.g. iRoddier et al.ll 19961) . If the 
system is edge-on, the disk itself can be seen as a dark lane 
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on the equatorial plane (e.g. HH 30: iBurrows et alJ (119961) . GM 
Auriga: IStapelfeldt & The WFPC2 Science Teaml l ll997l) ). 

Most of pre-main-sequence stars also show an ejection 
of matter. In several sources, optical linear jets are observed 
through the emission of fo rbidden atom i c line s. Interestingly, 
iBacciotti etafl (120031) and ICoffev et all ( 12004 have detected 
possible rotation signatures in at least four jets (DGTau, TH28, 
RWAur, LkHa321): They have observed systematic velocity 
differences between opposite sides of the jet axis. Indeed, ro- 
tation of the jet is predicted by magnetocentrifugal model s 
of accretion-ejection ( iKonigl & Pudritzil2000t IShu et al 1 120001 
through the coupling of the outflowing gas to the helicoidal 
structure of the magnetic field. In younger, more embedded 
sources, large molecular outflows are observed, mainly through 
the em ission of the rot ational lines of CO and its isotopo- 
logues (lArce et al.ll2006l and references therein). Those molec- 
ular outflows are usually identified with ambient molecular gas 
that has been entrained and thus put into motion by the under- 
lying protostellar jet. In any case, it is particularly interesting 
to check whether a rotation component is present in molecular 
outflows. 

In this context, one of the most interesting objects is the 
young PMS star HH 30, located in the Taurus molecular cloud 
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Table 1. Observation parameters. 
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a The rms noises are computed from the channel maps featured in Fig.|2|and|3| i.e. using natural weighting and the Sv channel resolutions of the 

maps. The original correlator resolution are between 3 and 6 times higher than the one shown here. 

b This is the total on—source observing time of useful data as if all observations were done with 5 antennas. 



at a distance of ~ 140 pc. iBurrows et alJ ill 9961) obtained spec- 
tacular images with the Hubble Space Telescope that revealed a 
flared edge-on disk of diameter ~ 450 AU, separated in two parts 
by a dust lane, and a highly-collimated jet emanating out of the 
central region. Since then ; the HH 30 disk and jet have been stud- 
ied in many details (e.g. Bacciotti et alJl!999t IStanelfeldt et all 
19991 IWood et al J 120001 ICotera et alJl200U IWpod et alJl2002t 
Watson & Staoelfeldtl 120041) . IStanelfeldt et alJ dl999l) reported 
variability and asymmetry in the emission from the circumstel- 
lar disk by comparing HST optical observations performed be- 
tween 1994 and 1998; the photometric variability of the star and 
its pos sible effects on the disk was investigated by I Wood et all 
J2000I This source is a unique candidate to probe the complex 
interaction between the disk, the jet and the outflow. 

However, HH 30 has been poorly investigated at mm wave- 
lengths so far. This is probably due to the confusion arising from 
the HH30 parent molecular cloud, which increases the com- 
plexity of any detailed observations using the emission of CO 
isotopologues. The 12 CO emission of the cloud is expected to 
be optically thick, hence hiding significantly the disk emission. 
BIMA observations of the 13 CO J=l-0 emiss ion in the environ- 
ment of HL Tau (Fig. 2 of IWelch et ail200f3l shows that HH 30 
is located in the south edge of a shell driven by XZ Tau. Still, 
the remarkable properties of HH 30 makes it a textbook case for 
the study of the disk/jet system in young PMS stars, and we have 
therefore decided to use the IRAM Plateau de Bure interferome- 
ter to perform high angular resolution (~ 1") observations of the 
emission of the circumstellar material surrounding HH 30 in sev- 
eral CO isotopologues lines. In this paper we present the results 
of this study and the analysis of the thermal dust and molecular 
line emissions of the disk and the outflow. 

2. Observations and data reduction 

2.1. 12 CO, HCO + and continuum at 3.4 and 1.3 mm 

First PdBI observations dedicated to this project were carried out 
with 5 antennas in BCD configurations (baseline lengths from 
24 to 280 m) during the winter 1997/1998. The observations 
were performed simultaneously at 89.2 GHz (HCO + J=l-0) and 
230.5 GHz ( 12 CO J=2-l). One correlator band of 10 MHz was 
centered on the HCO + J=l-0 line. Another band of 20 MHz was 
centered on the 12 CO J=2-l line. Finally, two bands of 160 MHz 
were used for the 1.3 mm and 3.4 mm continuum, respectively. 

Those observations were part of the surve y of pr otoplane- 
tary disks which is described by ISimon et alJ d200(Jl . The in- 
terferometer was operating in track sharing mode, i.e. each 8h- 



track was shared between several close-by sources. The total 
on-source observing time dedicated to HH 30 was about 6 hrs. 
The rms phase noises were between 8 and 25° at 3.4 mm and 
between 15 and 50° at 1.3 mm. This introduced position errors 
of < 0.1". The seeing, estimated from observations of the cali- 
brators, was ~ 0.3". Typical resolutions are 3.6" at 3.4 mm and 
1.4" at 1.3 mm. 

2.2. 13 CO, C ls O and continuum at 2.7 and 1.35 mm 

As a follow-up, we carried out observations of 13 CO and C 18 
at PdBI with 5 antennas in CD configuration (baseline lengths 
from 24 to 176 m) in December 2000, February 2001 and March 
2002. The C ls O J=l-0, 13 CO J=l-0 and 13 CO J=2-l lines 
were observed simultaneously using the 3 mm receiver (tuned at 
109.9 GHz) and the 1 mm receivers (tuned at 220.4 GHz). Three 
20 MHz correlator bands were centered on the C 18 J=l-0, 
13 CO J= 1-0 and 13 CO J=2-l lines. Two bands of 320 MHz were 
used for the 1.35 mm and 2.7 mm continuum, respectively. 

The rms phase noises were between 10 and 25° at 2.7 mm 
and between 20 and 50° at 1 .35 mm. The total on-source observ- 
ing time was about 14 hrs. Typical resolutions are 3" at 2.7 mm 
and 1.7" at 1.35 mm. 

2.3. Data reduction 

All dat a were reduc ed using the gildas 1 softwares supported at 
IPsAM dPetvl 120051) . Standard calibration methods using close 
calibrators were applied to all the PdBI data. Images were pro- 
duced using natural weighting of the visibilities for the line maps 
and robust weighting for the continuum maps. 

The 1.30 mm and 1.35 mm continuum observations were 
merged (after flux correction using a spectral index of 2.22, see 
Section l3~3l to produce a map of higher signal-to-noise ratio in 
Fig-H The channel maps presented in Figs.[2]and[3]include the 
continuum emission because it is in this case weak compared 
to the line emission. It has nonetheless been subtracted before 
fitting the disk (cf. section l4~Tl . 



3. Results 

All observed lines were detected although with very different 
strengths. 12 CO J=2-l, 13 CO J=l-0 and 13 CO J=2-l emissions 

1 See http://www.iram.fr/IRAMFR/GILDAS for more informa- 
tion about the gildas softwares. 
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Fig. 1. Left panel: C ontours of the CO J=2-l emission of HH 30 are plotted over a composite HST image (670 and 787.7 nm, see 
iBurrows et al Emission is only integrated at extreme velocities (< 4km s 1 and > 1 1 km s ') to avoid contamination by the 

disk emission. Contours are at the 3-sigma level (48 mJy/Beam) and the spatial resolution is 1.23 x 0.75" PA 31° (compared to a 
pixel size of 0.1" for the HST image). Medium panel : Superimposition of 3-sigma level contours 1) of the outflow emission as 
traced by the 12 CO(2-l) emission integrated at extreme velocities (cf. left panel) and 2) of the blue-shifted from 4.6 to 7.2 km s 
(blue dotted line, 29 mJy/Beam) and red-shifted from 7.2 to 9.8 km s _1 (red dashed line, 29 mJy/Beam) emission of the disk as 
traced by the 13 CO(2-l) line. The spatial resolution of the 13 CO map is 1.57 x 1.35" PA 25°. Right panel: Superimposition of the 
3-sigma level (1.4 mJy/Beam) of the merged 1.30 mm and 1.35 mm continuum emission over the same composite HST image. The 
spatial resolution is 1.29 x 0.91" PA 33°. 

Table 2. Continuum fluxes and sizes as a function of the frequency. 
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are strong while only weak C O J=l-0 and HCO + J=l-0 emis- 
sions are observed. Continuum emission at 3.4 mm, 2.7 mm, 
1.35 mm and 1.30 mm is detected, the source being resolved 
only at 1.30 mm. Fig. n summarize our molecular and contin- 
uum observa tions of HH 30 super imposed over the well-known 
HST image dBurrows et alj[l996l) . This figure clearly suggests 
that 1) the outflow is detected only through the 12 CO emission, 
2) the 13 CO J=2-l line mainly traces the rotating disk (the ve- 
locity gradient along the major axis of the optical disk is a direct 
signature of rotation), and 3) the mm continuum emission, cen- 
tered on the optical dark lane, is thermal and originates from 
the disk. Fig.[2]and[3]present the channels maps of the four ob- 
served CO lines. Those channel maps are centered around the 
HH 30 systemic velocity of the disk: 7.25 km s (cf. Section^. 
Only channels where signal is detected are shown, i.e. from -0.2 
to 13.8 km s" 1 for 12 CO J=2-l and from 4.9 to 10.3 km s" 1 for 
the other CO isotopologues. As a reference, the optical dark line 
and jet directions as seen on the HST image are sketched on 
those figures with orthogonal lines. 



3. 1 . The molecular cloud 

13 CO J=l-0 channel maps (Fig. |2]b) show emission centered 
on the optical disk on top of large structures, unrelated to the 
outflow and/or the disk. Those large structures are clearly seen 
on the 13 CO J=l-0 emission in the same velocity interval (i.e. 

6 to 8.4kms _1 ), both in the high spatial resolution PdBI data 
(Fig. |5]b) and in medium s patial resolution bu t wide field-of- 
view BIMA data (Fig. 2 of IWelch etaT1l2000t) . However, they 
are not detected on the I3 CO J=2-l channel maps (Fig.[5]a). We 
thus identify those large structures with the surrounding molec- 
ular cloud. As the disk is expected to be hotter than the molec- 
ular cloud, the 13 CO J=2-l line is mainly originating from the 
disk while the 13 CO J=l-0 emission shows both contributions. 
Moreover, the molecular cloud 12 CO J=2-l emission is ex- 
pected to be optically thick, hence hiding any emission from the 
disk/outflow system. Indeed, I2 CO J=2-l channels from 6.2 to 

7 km s _1 are devoid of signal in the disk/jet directions on Fig. [3] 
The interferometer has probably filtered out this emission be- 
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Fig. 2. Line emission of the CO isotopologues in HH 30. From top to bottom: 13 CO J=2-l, 13 CO J=l-0, and C ls O J=l-0 emissions. 
The channel maps are centered near 7.25 km s , the HH 30 systemic velocity. Blue and red contours respectively indicates the blue 
and red-shifted channels. Channel width is 0.32 km s~'. Plain and dotted lines respectively show positive and negative contours. 
Contour spacing corresponds to 3 <x for 13 CO J=l-0 and J=2-l and 2 cr for C 18 J=l-0. Noise levels and spatial resolutions may 
be found in Tabled The cross indicates the position and orientation 1) of the continuum emission at 1.3 mm and, 2) of the optical 
jet. 
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Fig. 3. Same as Fig.|5]for the emission of the 12 CO J=2— 1 line. Spectral resolution is 0.4kms _1 and contour spacing correspond to 
3 <t (cf. Table ^for more details). 



cause the size scale of the optically thick molecular cloud emis- 
sion is larger than typically half the interferometer primary beam 
(the well-known short-spacings problem). Finally, the emission 
of C ls O J=l-0 (the most optically thin line of this study) is very 
weak and peaks from 6.2 to 6.8 km s _1 (Fig.|2]c). 

3.2. The molecular disk 

Contrary to 13 CO J=l-0 line, the 13 CO J=2-l emission shown in 
Figgis clearly dominated by the circumstellar disk. The disk is 
seen in the 13 CO J=2-l emission from ~5 to 10 km s _1 . The disk 
is also detected in the same velocity range in the 12 CO J=2-l 
channel maps (cf. Fig.[3}- 

Fig.|4]a shows a Position- Velocity plot along the disk axis for 
12 CO J=2-l and 13 CO J=2-l. We overplotted 1) the direction of 



star and the systemic velocity as two orthogonal blue lines and 
2) the curves of the theoretical Keplerian velocity for a 0.45 M 
star. Note that this is just an illustration: The stellar mass value 
and uncertainty are determined through the detailed disk analy- 
sis described in section l4~2l A clear signature of the Keplerian 
rotation of the disk is visible. However, some confusion exists 
from the parent cloud. For 12 CO J=2-l, this confusion is so im- 
portant that the disk emission stays undetected due to the molec- 
ular cloud between 6.2 and 7 km s _1 . For 13 CO J=2-l, the confu- 
sion in the disk direction is important only for the 6.83, 6.51 and 
6.19kms~ 1 channels and minor elsewhere. The detailed mod- 
elling of the HH 30 circumstellar disk described in section l4~2"l is 
thus mainly based on the 13 CO J=2-l data in the velocity range 
devoid of confusion by the parent cloud. 
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Fig. 4. Position-velocity diagrams along the disk axis of the 
12 CO J=2-l (left panel, contour spacings set to 90 mJy/Beam or 
~ 3cr) and 13 CO J=2-l (right panel, contour spacings set to 60 
mJy/Beam or ~ 3cr) emissions. The blue horizontal and vertical 
lines respectively indicate the systemic velocity and the center 
of the disk as seen in the continuum. The red curves show the 
theoretical Keplerian velocity for a 0.45 M Q star. 
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Fig. 5. Left: Integrated emission of the HCO + J=l-0 line in the 5 
to 8 km s _1 velocity interval. Contour spacing corresponds to 1 <x 
(15 mJy/beamkms -1 ). The cross indicates the position of the 
continuum emission at 1.3 mm. Right: HCO + J=l-0 spectrum 
of the flux integrated over the ellipse shown as a red dash contour 
on the left panel (size: 10" x 5", with a major axis at PA: -50°). 
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With a critical density of ; 
HCO + J=l-0 transition is expected to be thermalized in proto- 
planetary disks. This transition is also easily observed in molec- 
ular clouds with a sub-thermal excitation, and often detectable 
in molecular outflows. However, the HCO + J=l-0 emission in 
HH 30 only comes from the disk as illustrated by Fig. [5] 

3.3. The dust disk 

We obtained continuum maps at four different wavelengths: 

3.4, 2.7, 1.35, and 1.3 mm. All of them show emission coin- 
cident with the HH 30 protostar position. This strongly suggests 
we are observing thermal emission from the dusty disk. Only 
the merged 1.35 mm and 1.30 mm continuum map is shown 
in Fig. U] As a first analysis, we fitted in the uv plane a bi- 
dimensional Gaussian function in each data set. Results are dis- 
played in table [5] The continu um fluxes agree with previ ous 
measurements made at OVRO (IStapelfeldt & PadgetJll999h . A 
fit of a power law, S v = S ioo(v/100GHz) Q ', through the measured 
fluxes yields a = 2.22 + 0.08 and S m = 3.19 ± 0.18 mJy. Low 
frequency continuum emission is unresolved while a Gaussian fit 
through the 1.35 mm and 1.30 mm merged data set gives a size 
of 1.26 ± 0.1 1 x 0.58 + 0.09" with a major axis at PA -50 ± 5° 
(i.e. a disk axis at PA 40 + 5°). This orientation is consistent with 
the orientation of the dark lane in the optical images (disk axis 



at PA 32.2 ± 1.0°. lBurrows et alJl99rl . Assuming a distance of 
140 pc, the fitted linear continuum FWHM is ~ 175 x 80 AU. 

3.4. The molecular outflow 

The outflow is only detected in the 12 CO J=2-l emission. 
However, as said before, there is confusion with the disk and 
parent cloud emission at velocities close to the HH 30 systemic 
velocity. Only the channels at extreme velocities (< 4 km s _1 and 
> 11 kms -1 ) can be attributed to the outflow without ambigui- 
ties. The integrated image from the outflow and the disk are pre- 
sented in Fig.|6] showing that the outflow contributes about twice 
as much to the total flux as the detected disk and parent cloud. 

Fig. |3 clearly shows that the CO outflow is essentially one- 
sided. While the North-East lobe is prominent, there is no emis- 
sion from South- West, except at low level (3 to 6 sigma) in the 
velocity range 7.8 to 8.6 km s _1 . In this northern lobe, the 12 CO 
emission delineates a conical structure, whose apex coincides 
(within the resolution of these observations) with the star loca- 
tion. The semi-opening angle is ~30°. The highly-collimatedjet 
seen by the HST is located precisely on the axis of this coni- 
cal 12 CO outflow, in agreement with the usual disk/jet formation 
paradigm. 

As expected from the edge-on geometry of the disk, the out- 
flow lies almost perfectly in the plane of sky. This is indicated by 
the outflowing gas being observed at both blue and red-shifted 
velocities in the same (northern) lobe, as well as by the quite 
low maximal velocities (only ~6kms _1 ). Indeed the extreme 
velocities, 1 and 13.2 km s~', are almost exactly symmetric 
from the systemic velocity determined from 13 CO (7.25 km s ). 
Assuming the 12 CO emission arises from outflowing gas along 
the cone wall, with a constant velocity expansion, we can de- 
rive an outflow velocity of 12 ± 2 km s _1 and an inclination of 
^0 + 2°. 

Although the molecular outflow is definitely observed as a 
conical structure emanating from the inner part of the disk, CO 
"clumps" also exist in projection along the jet axis, e.g. 1" NW 
off the source as illustrated in Fig. [7] 

4. Disk Model 

4.1. Analysis method 

A standard disk model JPringlelll981h was used to quantify the 
properties of the HH 30 circumstellar disk. This model assumes 
an axi- symmetrical geometry, local thermodynamic equilibrium 
and the absence of a vertical gradient of temperature. Turbulence 
is phenomenologically introduced by adding a turbulent width 
Av to the thermal width. Radial evolutions of temperature, sur- 
face density and velocity are assumed to be power law normal- 
ized at 100 AU: T = Tm (r/100AU) -? , 2 = 2 100 (r/100 A\JT P 
and V = Vioo (r/100 AU)" V . A value of 0.5 for the velocity scal- 
ing exponent v implies that the disk rotation is Keplerian. We 
further assume that the density profile follows a Gaussian distri- 
bution of the height z above the disk plane with the scale height 
H(r) being a power law of r 



n(r, z) — n(r) exp < - 



H(r) 



and n(r) = 



H = H m (r/100 ALT and n = n m (r/lOOAU) 

which implies i 
imply h = 1 + \ 



which implies s — p + h ( hydrostatic equilibri um would further 
q/2, see lDartois et al.l2003l) . 
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Fig. 6. Integrated emissions of the 12 CO J=2-l line in three contiguous velocity intervals: 1) blue-shifted velocities where only the 
outflow emits (left panel), 2) near-systemic velocities where confusion arises between the outflow, the disk and the cloud (medium 
panel), 3) red-shifted velocities where only the outflow emits (right panel). Maps have been rotated by 30° clockwise around the 
disk center to bring the optical jet direction in the vertical direction as this eases the extraction of position-velocity diagrams. This 
defines a new coordinate system, named x, y hereafter. In addition to the horizontal-vertical cross featuring the dark line and the 
optical jet, a second cross is sketching a cone of 30° half-opening angle. The contour spacing corresponds to 3 cr (i.e. from left to 
right: 0.16, 0.18 and 0.16 Jy/Beamkms -1 ). Fluxes estimated inside polygons following the 2 <x contour (dashed lines) are from left 
to right: 6.8, 7.0 and 8.3 Jy km s _1 . 
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Fig. 7. Spatial and kinematical properties of the CO J=2-l 
"clump" observed along the jet axis: (a) Top left: xv position- 
velocity diagram averaged on 5" < 5y < 8". (b) bottom left: 
Integrated intensity in the [0,3] and [10.5,13.5] km s _1 velocity 
ranges, (c) Bottom right: vy position-velocity diagram averaged 
on -0.5" < 5x < 0.5". (d) Top right: Spectrum integrated in the 
-0.5" < Sx < 0.5" x 5" <5y< 8" rectangle. 

The model parameters are: D the distance from Earth, Vlsr 
the systemic velocity, PA the plane-of-sky orientation, R out the 



disk outer radius, Av the turbulent line width, and the parame- 
ters of the power laws ({Vioo,W, \T\oo,q}, {HiooJi} plus a pair 
[Sioo, p) for each molecule transition and for the dust). Since the 
whole model is described by a limited number of parameters, it 
is possible to perform a x 1 minimization in order to derive the 
best-fitted values of each parameter. The fitting is performed in 
the uv plane to avoid any error introduced by the imaging and 
deconvolution stages. The^ 2 is thus defined from the difference 
between the observed and predicted visibilities: 

2 

, v-i model(M, v, V) - observed(M, v, V) 

where py„ v i is the noise associated to each visibility (see 
iGuilloteau & Dufrevlll998T) . The minimization were done after 
subtraction of the continuum visib ilities from the line uv ta- 
ble (for details, see lPietu et alJ2006l) . 

The analysis results are summed up in Tables [3] and |4] The 
comparison between our best model and the observations is 
shown in Fig.|F]for the 13 CO J=2-l line. 

4.2. 13 CO lines 
4.2.1. Fit difficulties 

The near edge-on geometry of HH 30 makes this source a spe- 
cial case. First, the model must properly sample the disk thick- 
ness to adequately model the line emission. Second, the num- 
ber of independent data points is actually small since the emis- 
sion is unresolved perpendicular to the disk plane. This implies a 
significant degeneracy between parameters Tioo, q, 2, p, Av and 
the scale height H(r). In fact, the 13 CO J=2-l line intensity 
depends weakly on the temperature in the expected tempera- 
ture rang e for HH30, i.e. a round 20 K for the linear scales 
sampled foartois et al.ll2003l) . At any (non zero) projected ve- 
locity Vobsi the line flux will mostly come from a radius r = 
100(Vioo/i'obs) 2 , and the total emitted flux will be essentially 
proportional to l.(r)H(r)5V, where 6V is the local intrinsic line 
width. 

A further difficulty is the existence of contamination by the 
molecular cloud, and perhaps also by the outflow. We avoided 
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Fig. 8. Channel maps of the I3 CO J=2-l emission. Observations, 
best model and difference (= observation-model) are plotted 
from left to right. The coordinate system is native (i.e. equato- 
rial). The spatial resolution is 1.8x1 .6" at PA of 26°. All contours 
are equally spaced by a value of 45 mJy/Beam (2.5 cr at the spec- 
tral resolution of 0.64 km s~'). Blue and red contours respec- 
tively indicate the blue and red-shifted channels while the black 
contours indicate the systemic velocity channel at 7.25 km s~'. 
Plain and dotted lines respectively show positive and negative 
contours. The cross indicates the disk center as measured in 
the 1 .3 mm continuum emission. A geometrical asymmetry (5cr- 
contours in the difference maps) is detected in the blue-shifted 
part of the disk, close to the systemic velocity. 



contamination by the molecular cloud by ignoring the velocity 
range from 6.2 to 6.8 km s _1 in the determination of the x 1 - Any 
contamination by the outflow would separate from the disk only 
spatially, not in the velocity space. As the fit is done in the uv 
space, contamination by the outflow is unavoidable. Although 
the outflow contribution is most likely small, because emission 
at large velocities from the outflow region is not detected in the 
I3 CO J=2-l channel maps (see Fig. [HJ, it may slightly bias the 
inclination towards lower values. 



4.2.2. Fit steps and results 

We have performed x 1 global minimizations using the follow- 
ing set of parameters: PA, i, R out to constrain the geometry, and 
Visr, Vioo; v, to constrain the velocity power law, using fixed val- 
ues for the other parameters (see Table [5}. The systemic veloc- 
ity value of 7.25 km s~' quoted before comes from this step. We 
verified that the rotation is Keplerian on the pair of parameters 
(Vioo, v )- This again suggests that any remaining confusion due 
to the outflow or cloud is small. We found HH 30 to be essen- 
tially edge on, with a best fit inclination as 81°. We nevertheless 
fixed the inclination to 84° in all other fits to avoid the possible 
inclination bias due to the outflow contamination (see above). 
We verified that all other parameters are essentially independent 
on this small variation of the inclination. 

The degeneracy between the remaining parameters cannot 
be fully removed with the available data. We thus assumed the 
height law parameters to be Z/ioo = 22 AU, and h = 1.25. 
We then fitted the J=l-0 and J=2-l lines of 13 CO simulta- 
neously. We find for the temperature law: T = 12 ± 1 K , 
q = 0.55 ± 0.07. The associated 13 CO surface density at 100 



AU is about 9 10 cm (with a factor two uncertainty), but 
the surface density power law index p is less well constrained: 
p - 1.5 ± 1. 

Temperature and density power laws are less certain because 
1) the simultanesous fit is only valid in the absence of strong 
vertical temperature gradient and 2) the constraint uses only a 
few channels from the J=l-0 line of 13 CO where the contami- 
nation by the parent cloud is thought negligeable. We neverthe- 
less stress that our assumptions have essentially no influence on 
the main disk parameters, Vioo, *, and R out , which are well con- 
strained by the current data. 



4.3. HCO + J=1-0line 

HCO + emission is detected from the disk, but with a relatively 
low signal-to-noise ratio. There is no evidence of HCO + emis- 
sion either from the cloud or from the outflow. Accordingly, the 
HCO + data were analyzed in a very similar way to the 13 CO data, 
but all channels were included in the analysis. The results (e.g. 
outer radius) agree within one sigma with those obtained from 
I3 CO, but the error bars are large. In practice, only the HCO + 
column density is constrained from these observations. Using 
p = 1.5 and R oat = 400 AU, we find 2[HCO + ] = 7.2 ± 1.5 10 12 
crrT 



at 100 AU. 



4.4. Continuum Emission 

We have high resolution continuum observations at four differ- 
ent frequencies. As the disk dust thermal emission is at least 
marginally resolved and partially optically thin, we may obtain 
an independent estimate of the outer radius and the surface den- 
sity. Matching the surface density found in this process with the 
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one derived from the 13 CO data enables to derive the abundance 
of 13 CO. 

We used the disk model derived from the 13 CO data and the 
four continuum measurements to constrain the dust emissivity 
index, /3, and the surface density Eioo assuming the emissivity 
law is given by k(v) = O.lCv/lO^Hzy 3 cm 2 g _1 . We made a global 
X 1 fit using the following four parameters: R oat , fl, Zioo and p. 
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Table 3. Comparison of mm and optical/NIR properties. The upper part of the table shows parameters whose value is robustly 
deduced from the millimeter data (geometry, velocity law and stellar mass) while the bottom part displays parameters less well 
constrained due to the limited vertical resolution of this edge-on disk. Most parameters are deduced from the analysis of the 13 CO 
emission, except for the parameters of the H2 surface density law deduced from the continuum emission. 





This work ( 13 CO & continuum data) 


IBurrows et alj 
119961) HST 


ICotera et alj 
<20011) HST 


Wood et al. 

<2002t) 
HST+SED 


Assumed Distance 


D (pc) = 140 


140 


140 


140 


Systemic velocity 
Orientation (disk axis) 
Inclination (disk plane) 


V, sr (km.s- 1 ) = 7.25 ±0.04 
PA O = 32 ±2 

(81 ±3 (Best fit) 
' (84 ±3 (Canonical fit) 


32 
82.5 


32 

84 (assumed) 


32 

84 (assumed) 


Outer radius 
Turbulent linewidth 


flout (AU) = 420 ±25 
Av (km.s- 1 ) = 0.23 ±0.03 


(250 (Canonical fit) 
\425 (Best fit) 


200 (assumed) 


200 (assumed) 


Velocity 1 

Velocity at 100 AU 
Velocity exponent 
Stellar mass 


aw: V(r) = Vino I — 1 

U00AU/ 

Vioo (km.s- 1 ) = 2.00 ±0.09 
v = 0.50 ±0.06 
M»(M )= 0.45 ±0.04 


0.67 (assumed) 


0.5 (assumed) 


0.5 (assumed) 


Scale Heig 

Scale Height at 100 AU 
Scale Height exponent 


/ r \- h 

it law: H(r) = H lm 

U00AU/ 

H 100 (AU) 22 

h 1.25 (assumed) 


22 (= 15.5 V2) 
1.45 


21 (= 15 V2) 
1.29 (= 58/45) 


24 (= 17 V2) 
1.25 (assumed) 


Temperatur 

Temperature at 100 AU 
Temperature exponent 


elaw: T(r) =T l0 o( mM] ) 

* 100 UV) ±i 

q = 0.55 ±0.07 


34 (Disk surface) 
0.1 (= 3-2/;) 


0.4 (= 3 - 2h) 


0.5 (= 3 - 2h) 


H? Surface De 
Surf. dens, at 100 AU 

Surf. dens, exponent 
Density exponent 


.nsitylaw: S(r) = E 100 ( ^ A J 

E100 (cm" 2 ) = 3.6 ±0.6 10 22 
E100 (g-cnr 2 ) = 0.16 ±0.03 
1.0 

S =a 2.2 


5.5 10 21 
0.024 
0.75 (= s - h) 
2.2 


1.08 (= s-h) 
2.37 (assumed) 


1 (= s - h) 
2.25 (assumed) 



Table 4. HH 30 mm dust properties. The second column displays the best fit results while the third column shows the canonical fit 
results {e.g. obtained with a fixed outer radius of 420 AU, see text for details). The total mass is calculated assuming a gas-to-dust 
ratio of 100 and using the kinetic temperature derived from the 13 CO analysis (Tioo = 12 K, q = 0.55). 





Dust: k v = 


K ° X ^10 12 Hz^ 




Absorption law 


k (cm 2 .g-') = 


0. 1 (assumed) 




Dust exponent 


P = 


0.4 ±0.1 


0.5 ±0.1 


Absorption at 230 GHz 


k(230) (cm 2 . g- 1 ) = 


0.055 


0.048 


Dust disk size 


Rj (AU) = 


145 ±20 


420 (assumed) 


Surface Density (H2) 


Z100 (cnT 2 ) = 


8.6 10 22 ±1.410 22 


3.6 10 22 ±0.6 10 22 


Exponent 


P = 


±0.5 


1.0 ±0.1 


Total mass 


M ilA (M ) =; 


2.7 10~ 3 ±0.4 10" 3 


4.8 10~ 3 ±0.8 10~ 3 
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We assumed identical gas and dust temperature at 100 AU. 
The best fit is obtained for a uniform disk (p = 0) of small radius 
(150 + 20 AU) and total mass 2.7 10~ 3 M , but solutions with a 
decreasing surface density and larger radii remain acceptable. 
Assuming the same outer radius than for CO leads to p = 1.1, 
and a surface density at 100 AU of 3.6 + 0.3 10 22 ctrT 2 , which 
corresponds to a mass of 4.8 10~ 3 M Q . Note that the mass scales 
inversely with the assumed dust temperature. In all cases, we 
frndj3 = 0.4±0.1. 

5. Outflow Model 

5. 1 . Model description 

While a geometrically thin disk with power law distributions of 
the radial dependence of density and temperature constitute a 
good description of circumstellar disks, no such paradigm exist 
for outflows. Outflow models still rely on ad hoc parametriza- 
tion, with basic ingredients often differing from source to source. 

The 12 CO outflow of HH 30 presents a remarkably simple 
structure, both in terms of morphology and velocity distribution 
that we modelled as sketched in Fig.|9] 

Geometry - We assume a perfect conical geometry, with an 
semi-opening angle £> max . The cone axis is inclined by an angle i 
to the line-of-sight. 

Velocity distribution - The velocity of the outflowing gas 
is assumed to have two components at any position: a radial 
component, i.e. this component vector is pointing outward of the 
central source position and an azimuthal component (to trace a 
possible outflow rotation), i.e. a vector component tangent to the 
circle defined by the intersection of the outflow cone with a plane 
perpendicular to the jet axis. 

Density distribution and Excitation - Our model is assum- 
ing optically thin isothermal emission: The emissivity of the gas 
is proportional to the local density and the brightness distribution 
is proportional to the column densities. The model is computed 
in arbitrary units, and the peak of the resulting image is further 
scaled to the peak of the observed maps. Our approach is thus 
to focus on the kinematics and morphology of the outflow rather 
than on the excitation conditions. We assumed that the density 
is located mainly along the cone edges with a constant Gaussian 
width scale Wq and a Gaussian decrease of the density with the 
distance to the star along the jet axis (Ho being the height scale). 

The model is computed on a sufficiently fine grid, and is fur- 
ther convolved with the clean beam of the observations, in order 
to mimic the same angular resolution. 

5.2. Best model 

Fig.llOlshows a comparison of the 12 CO data and of our "best" 
model. This comparison is composed of images of the integrated 
emission and of two position-velocity diagrams (perpendicular 
to and along the jet axis) and of the spectra integrated over the 
outflow. Complete sets of channel maps and position-velocity 
diagrams are presented in Fig. ^2 D3 an d El (electronic ver- 
sion only). While the 12 CO emission clearly reveals the outflow 
structure, it also includes a contribution from the rotating disk 
and the parent cloud. Hence, we blanked out the velocities be- 
tween 4.6 and 9.8 km s _1 in the data cube in order to make mean- 
ingful comparisons between data and model. While this is clear 
on the position-velocity diagrams and the integrated spectrum 
of Fig. El note mat me integrated emission map shown as the 
top, left panel also excludes this velocity range. Fig. I14landll5l 
(electronic version only) are variations of Fig.El wn ich illustrate 




Fig. 9. Sketch of our outflow model. The emitting gas is con- 
fined in a layer near the edge of a conical structure (grey zone). 
The dotted vertical line indicates the plane of the sky. The south- 
ern part of the disk is pointed t oward us in agreeme nt with the 
scattered light observations (e. e. iBurrows et al.ll996l) . 



how the different parameters of the outflow model influence the 
observables. 

We did not develop any x 1 fitting procedure. We took in- 
stead the following steps to successively find a plausible range 
of parameters representing the data. The systemic velocity was 
assumed to be the same for the outflow as for the disk. The po- 
sition angle was tuned to give the best averaged horizontal sym- 
metry of the outflow on each side of the jet on the channel maps 
of Fig.[^ The opening angle of the cone was measured on the 
vy position-velocity diagrams (Fig.ll3> as the averaged slope of 
the minimum y value where emission is detected as a function of 
the x plane. The magnitude of the turbulent line width was de- 
duced from the slope of the outer wings of the integrated spec- 
trum. Indeed, Fig. 1151 shows that the absence of turbulent line 
width makes infinitely sharp outer wings. Once the systemic ve- 
locity is fixed, only the inclination on the plane-of-sky ensure 
the correct velocity centering of the emission features in the vy 
and xv position-velocity diagrams. Once the cone opening an- 
gle and the turbulent width are fixed, only the radial velocity 
can ensure the correct velocity width in the vy and xv position- 
velocity diagrams. The density width and height scales were fi- 
nally tuned to reproduce the width and height of the horn in the 
vy position-velocity diagram. Table[5]lists the parameters of the 
"best" model and associated possible ranges of parameters. We 
checked that the vertical limitation of the outflow is mainly phys- 
ical, i.e. not produced by the short-spacing filtering of the inter- 
ferometer. Finally, rotation was searched for but not found (see 
next section). 

The model depicted in Figs.ll0ltoll3land Table [5] nicely re- 
produces the main characteristics of the HH 30 outflow, both in 
terms of velocity and morphology. Hence, we conclude that, to 
first order, the outflowing gas in HH 30 forms a cone with a 
constant radial velocity distribution and no detectable ro- 
tation (see next section). The inclination on the plane-of-sky 
(constrained to be essentially -1° ± 1°) indicates that this cone 
is nearly, but not exactly, perpendicular to the disk. 

Nevertheless, a number of features are not reproduced by 
the model. First, our fully axisymmetric model can not repro- 
duce the various asymmetries seen: 1) Fig. 1131 illustrates that 
the cone opening angle is correct for negative 5x offsets but too 
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Table 5. Outflow properties. The values in this table are model 
dependent. Furthermore they do not results from a fit. They thus 
are very different in nature from the values in the tables describ- 
ing the disk which results from a validated fitting method. 



Observations 



Model 





"Best" Possible range 


Assumed Distance 


D (pc) = 140 (assumed) 


Systemic velocity 
Orientation 
Inclination 
Half-opening angle 


Via (km.s -1 ) = 7.25 ±0.04 
PA O = 32 ±2 
i(°)= -1 ±1 
e raax (°)= 30 ±2 


Turbulent linewidth 


Av(kms-')= 0.40 ±0.15 


Velocity law: 

Radial velocity 

Rotation velocity at 200 AU 


V r = v iad , V g = 0, V = V rot200 r* 

V rad (kms" 1 )= 11.5 ±0.5 
V ro , 200 (kms- | )= <1 


Volume density law: n c 
Width 

Scale Height 


<4(^^)H-(a)"} 

w (") = 0.3 < 0.3 
H (")= 3 ±0.5 



small for positive 5x offsets; 2) The position- velocity diagrams 
and the integrated spectra of the data clearly shows a brightness 
asymmetry between the [0.6,4.2] and the [9.8, 13.4] kms -1 ve- 
locity ranges. Second, the brightness distribution is more cen- 
trally peaked in the data than in the model (maybe due to hetero- 
geneous excitation conditions). Third, there is the high velocity 
"clump" 1000 AU from the star along the flow axis. The associ- 
ated emission shows a continuity both in the red and blue-shifted 
velocities with the cone emission that surrounds the optical jet 
(See the v- Sy diagram in Fig.0. This emission peak thus looks 
like a "clump" in the cone of outflowing gas, in contrast with 
standard CO "bullets" which show a distinct emission peak at 
vel ocities s ignificantly higher than the bulk of the outflow (See 
e.g. lBachiller et alll9 90)r 

Furthermore, because of the low inclination, our observa- 
tions may be insensitive to the component of motion parallel to 
the flow axis, because such a gas would appear at the systemic 
velocity and be indistinguishable from the surrounding cloud. 
Since we do not know the cloud extent, we can neither prove nor 
disprove that such a confusion is happening. 



5.3. Rotation? 

There is no strong prescription on the dependence of the rota- 
tion velocity with distance from the star and the jet axis. To 
test the possibility of rotation, we have introduced several ro- 
tation laws in the model: 1) constant velocity, 2) solid rotation 
V rot oc r, 3) Keplerian rotation, and 4) vortex rotation V rot oc 1/r 
(corresponding to angular momentum conservation in an ex- 
panding thin ring that would generate the outflow cone). In all 
those cases, rotation would manifest as a tilt of the ellipse of 
the (V, Sx) position-velocity diagram (second row of Fig. [TO] see 
also Fig. left column). We selected as upper limit to V rot the 
value that tilted the ellipse by 5°. In all cases, we can firmly rule 
out values of V rot larger than 1 km s _1 at 200 AU from the jet axis. 
In addition, Keplerian and vortex rotation would give rise to low 
level but detectable wings outside the displayed velocity range. 




5 10 
Velocity (Hei.tf" ) 



O 5 10 

Velocity (fcm.a ) 



Fig. 10. Comparison between 12 CO J=2-l observations (left col- 
umn) and our best model (right column) of the outflow. The 
velocity channels where confusion from the disk or the parent 
cloud exist have been flagged in the I2 CO J=2-l data cube. The 
integrated emission map, the position-velocity diagrams perpen- 
dicular and parallel to the jet axis and the spectrum integrated 
over the map are shown from top to bottom. The plots are shown 
in the rotated (Sx, Sy) coordinate system. The horizontal and ver- 
tical lines on the integrated emission map show the position of 
the cuts used to form the position-velocity diagrams (respec- 
tively Sx = -0.1" and Sy = 3.1"). The two horizontal or verti- 
cal lines in the position-velocity diagrams indicate the velocity 
range over which the emission have been integrated to form the 
top map (0.6 < V < 13.4 km s -1 ). 



Such wings are not detected in the complete data set, probably 
ruling out these kinds of rotation. 
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6. Discussion of stellar and disk parameters 

6.1. Stellar mass estimate, luminosity and spectral type 

Using a distance of 140 pc and the magnitude of the Keplerian 
rotation, we deduce a mass of 0.45 M for the young star. To use 
this stellar mass value in a distance independent evolutionary 
track diagram |log(L»/Mj), log(r eff )}, we also need the stellar 
luminosity value. The current estimates for HH 30 are uncertain 
beca use most of the stella r light is intercepted by the edge-on 
disk. iBurrows et alJ i 1996ft derived this stellar luminosity from 
scattered light emission at the disk surface, under the assumption 
of thermal equilibrium, and found (using our definition of H(r)): 

L = 1.0L (// 10 o/22A[/) 12 (M/0.67M ) 6 . 

Using flioo = 22 and M = 0.45 M Q , the corresponding luminos- 
ity is L — 0. 1 Lq, with at least a factor 2 uncertainty due to the as- 
sumed scale height . Based on the analy s is of the spectro-energy 
distribution (SED~). lKenvon et all i ll 9941) hav e estimated the stel- 
lar luminosity to be ~ 0.2 L Q . More recentlv. lCotera et alJll2001l) 
have estimated the luminosity to be between 0.2 and 0.9 L Q . 

Assuming L» = 0.2 L as the more probable value and using 
iBaraffe et alfi 19981) tracks (see Fig. 3 of lSimon et all2000l) . we 
find that the spectral type seems to be M3 with an age around 
4-5 Million years. Note however that the tracks are almost par- 
allel to the ordinate, i.e. log(L*/M 2 ), in this range of mass and 
luminosity. A relatively small error on the luminosity leads to 
an important uncertainty on the age. If we take the other end 
of the luminosity range (i.e. 0.9 L ), we get the same spectral 
typ e for an object wh ich is less than 1 Million year old. Using 
the lSiess etallfeOOOl) tracks, the spectral type would be M1-M2 
with an age between 4 and 1 Myrs for a lu minosity equal to .2 
an d 0.9 L Q , respective l y (see again Fig. 3 of lSimon et all20 00). 

lAppenzeller et alJ ( l2005h derive a spectral type of K7 for 
HH30, but caution that this re sults is based on a l imited wave- 
length coverage. From Fig. 3 o flSimon et allfeOOOl) . such a spec- 
tral type is incom patible with a mass of 0.45 M for the evo - 
lutio nary tracks of IBaraffe et all ( 119981) . IPalla & Stahleil ( I1999I) . 
and ISiess et all (|2000) which fi t all other sources. A K7 spec- 
tral type is compatible for the iD'Antona & Mazzitellil i 19971) 
tracks, but these fail to g ive agreement for the other stars. The 
lAppenzeller et alJ d2005t) result nevertheless suggest that a spec- 
tral type of M3 is unlikely. 

Adding the existence of a relatively powerful jet, it then 
seems reasonable to classify HH 30 as a young TTauri of stellar 
mass 0.45 M Q , spectral type ranging around Ml and age around 
4 Myrs. 

6.2. Disk Geometry 

Both the mm and the optical observations reveal a disk close to 
edge-on. Table [5] quantitatively compares the results from both 
wavelength ranges. Optical & NIR observations indicate that the 
disk is tilted along the l ine of sight with i ranging from 82 to 
84° llBurrows et alJ 19961 and references therein). In all cases, the 
PA is found to be 32° counted anti-clockwise from North. Our 
results agree well with these values. 

T he determina t ion of the outer radius i s more model depen- 
dent. ICotera et all ( 1200 ll) andlWood et al.1 d2002l) ha ve assumed 
R out = 200 AU while IWatson & Stapelfeldtl ll2004l) have taken 
R out = 250 AU. All these are consistent with the minimum ra- 
dius derived from the 1.3 mm continuum emissio n. However, in 
their best simulations (i.e. not the canonical one). lBurrows et all 



ill 9961) have found R oui ^ 440 AU, in excellent agreement with 
the value derived here from 13 CO J=2— 1. 

6.3. Disk Mass and Dust Properties 

Disk masses were previously derived by several authors from 
NIR and optical scattered light images l lBurrows et al.l fl996l 
ICoteraetal.ll2001l). som etimes combined with modelling of the 
SEP llWood et al.l2002l and references therein). 

Our derived disk mass depends on three assumptions: the 
temperature law, the surface density exponent p and the dust 
emissivity. In the optical and NIR, only the upper layers of the 
disk surface are seen. Therefore, a detail modelling of the surface 
density distribution cannot be achieved without some a priori 
assumptions even though th e width of the dark lan e can be used 
to estimate the opacity law. IBurrows et alJ i ll 9961) extrapolated 
a total disk mass of 6 10~ 3 M o , including a correction factor of 
1 5 to account for the likely existence of a vertical temperature 
gra dient. 

IWood et al.1 l l2002l) derive a mass of 1 .5 10" 3 M G . To do this, 
they use a SED fitting procedure and more sophisticated dust 
properties (e.g. with a grain size d istribution) leading to an emis- 
sivity of ~ 40cm 2 g 1 at 0.6 ^m. IWood et alJ ll2002l) dust prop- 
erty implies ft = 0.75 and an emissivity of 0.08 cm 2 g~' at 
1 .3 mm. The difference with our determination is thus mainly 
due to the assumed dust properties (although the temperature 
and disk size also have an influence). 

The value /3 = 0.4 - 0.5 is at the low end of the range en- 
countered i n proto-planetary disks, and very su ggestive of grain 
growth (e.g. lTesti et alJl2003t iNatta et alJl2004l) . As we have re- 
solved the disk in its radial direction, and since our model in- 
cludes the vertical distribution of the dust, this value of f3 is 
not biased by the contribution of an optically thick region 1 . 
However, the determination of ft depends on the assumption that 
the 3.4 mm flux is due to dust emission only. A small contam- 
ination by some free-free emission from the jet would drive ft 
towards lower values. Additional low frequency measurements 
are needed to evaluate this problem and confirm the low f3 value. 
iDraind ll2006i) demonstrated that, provided the maximum grain 
size exceeds about 3A (i.e. 1 cm in our case), such low /3 values 
can be obtained for power law grains size distribution n(a) oc cT y 
with exponent y ^3.3, since the apparent /3 index is linked to the 
material emissivity /3 e by /3 = (y - 3)/3 € . 

6.4. Molecular Abundances and Disk Temperature 

Matching all the surface density measurements indicate a 13 CO 
abundance of ^ 2 10~ 6 , and an HCO + abundance of 2 10~ 10 , 
with likely variations by a factor 2 between 100 and 400 
AU, since the column density law exponents do not necessarily 
match. An additional uncertainty of a factor 2 should be added 
due to the uncertainty in the dust emissivity at 1.3 mm. 

The temperature of 12 K derived from 13 CO data is signif- 
icantly lower than the value of 34 K derived by IBurrows et alJ 
(119961) . This is expected since the optically thin 13 CO emis- 
sion mostly traces the disk plane, while the scattered light traces 
the disk atmosphere. A decreasing temperature towards the disk 
plane is expected in a disk which is optically thick to the stel- 
lar radiation, but thin to its own re-emission. A rather surprising 
fact is that, despite this low temperature, the 13 CO abundance 

2 The disk is not resolved in height, however. If the dust was confined 
to a very thin layer (scale hight smaller than a few AU), and seen edge- 
on (i > 87°), an additional correction would be necessary. 
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suggests a low depletio n of CO compared to the surrounding 
Taurus molecular cloud dCernicharo & Guelinl fT987h. in contrast 
with what is found in most disks (e.g. lDartois et al.l2003l for DM 
Tau). 



7. Discussion of outflow parameters 

7.1. Outflow mass and energetics 

In the northern lobe of the outflow, the total 12 CO J=2-l line 
flux is 15.1 Jykms" 1 (using only velocities in the [0.6,4.2] and 
[9.8, 13.4] kms~' range and integrating spatially in the area de- 
fined in Fig. |6j. The 13 CO J =2-1 line flux from the same re- 
gion and velocity range is » 0.2 ± 0.05 Jy km s _1 . This implies a 
12 CO/ 13 CO ratio of -80. Thus the 12 CO emission is essentially 
optically thin. Further assuming it is at LTE allows us to de- 
rive the outflowing mass. A minimum flow mass of 2 10~ 5 M is 
found for temperatures in the range from 15 to 25 K, using a CO 
abundance of 10 4 . The outflow mass scales approximately as 
the temperature for larger temperatures. The volume density in 
the outflow can be derived from its total mass and an estimate of 
the volume it fills. Using the conical geometry described in the 
previous sections and a cone thickness of 45 AU (0.3") we ob- 
tain a mean density n(H2) — 8 10 4 cm~ 3 , high enough to validate 
the LTE assumption. 

The HH30 CO outflow presents an obvious, extremely 
strong asymmetry, since almost no signal is detected in its south- 
ern lobe. We estimate here the 12 CO J=2-l line flux to be less 
than about 1.5 Jykms 1 (over the same velocity range as that 
used in the northern lobe and in a symmetric region). This im- 
plies that this side of the flow is more than 10 times less mas- 
sive than the northern one. Note that, while the southern optical 
jet brightness is also weaker than its northern c ounterpart, the 
densi ty ratio has been estimated to be only ~ 2 dBacciotti et alJ 
H999T) . 

The physical parameters of the northern jet and outflow esti- 
mated in a region of ~ 6" or ~ 800 AU from the central source 
(corresponding to the size of the maps presented in this paper) 
are listed in Table [6] The values for the outflow were derived 
from the observations presented in this paper. Those for the jet 
have been estimated from the literature . The mean jet density 
was estimated to be around 10 4 cm' 3 bvlBacciotti et alJ ([1999). 
Using a jet radius of 0.09" dBacciotfi et alJI 19991) and a cylin- 
drical geometry 3 , the jet mass up to 800 AU is about 2 10~ 8 M Q . 
iBurrows et alJ i 19961) indicates that the jet knot velocity range 
between ~100 and 300 km s -1 , so we used a typical jet velocity 
of 200^8-'. 

7.2. Formation of the outflow - entrapment? 

The origin of the molecular outflows is a pending problem in 
star formation theories. On large scales (» 1000 AU), there 
is strong evidence that the observed outflows consist in ambi- 
ent molecular gas that has been put into motion by large bow 
shocks propagating down the underlying protostellar jets (e.g. 
iGueth & Guilloteaulll999l) . In this classical "prompt" entrain- 
ment scenario, the jet transfers momentum to the outflow via 
the bow-shock, and it is usually ass umed that the momentum 
flux is conserved in this process (e. g. iMasson & Cherninl 19921) . 
Interestingly, this is what is indicated by our first order estima- 
tion (see table[6}. 

3 The opening angle is only 2° beyond 100 AU, so the cylindrical 
approximation is appropriate 



Table 6. HH 30 jet and outflow parameters in the northern lobe, 
measured in a region of ~ 6" or ~ 800 AU from the central 
source (see text). 



Quantity 


Jet 


Outflow 


Mass ( Mq) 


210~ 8 


210~ 5 


Velocity (kms~') 


200 


12 


Dynamical timescale (yr) 


20 


320 


Mass flux (M G yr _1 ) 


1.0 10-' 


6.3 10- 8 


Momentum ( M Q km s _1 ) 


4.0 10~ 6 


2.4 10"* 


Momentum flux ( M Q km s~'/yr) 


2.6 10~ 7 


7.5 10~ 7 



Viscous entrainment alon g the jet edges naturally results in 
conical shape for the outflow lStahlerl[T994l) . However, the pre- 
dicted velocities are essentially parallel to the flow axis, while 
our observations are essentially sensitive to the radial motions 
because the flow axis lies in the plane of the sky. The compari- 
son is thus impossible in practice. The insensitivity of our obser- 
vations to the gas flowing near the direction of the jet axis also 
implies that the outflow momentum flux we derive above is only 
a lower limit. 

If the HH 30 outflow does indeed consist in accelerated am- 
bient gas, the lack of CO emission in the southern lobe would 
point towards different properties of the interstellar medium in 
this direction, that would strongly alter the outflow formation 
process. In this context, it may be worth mentioning two similar- 
ities between HH 30 and L 1 157 dGueth et aljll996l) . While the 
sampled scales are different (several 10 4 AU for L 1 157, only a 
few hundreds AU in the HH 30), both outflows exhibit 1) a more 
or less pronounced asymmetry between their lobes and 2) a con- 
ical geometry and a radial velocity near the launching region. 

7.3. Formation of the outflow - disk wind? 

One of the most striking results revealed by our data is the fact 
that the outflowing molecular gas is continuously collimated 
down to the very close vicinity of the star, at spatial scales that 
are smaller than the disk size (see Fig.^. This is somewhat diffi- 
cult to reconcile with the propagation of a large bow-shock (see 
previous section). But this is naturally explained if the observed 
CO emission arises from material that has been directly launched 
from the disk and evolves ballistically in the first few 100AU 
from the star. This hypothesis is also consistent with a coni- 
cal flow structure, with constant radial velocity that is observed. 
Recollimation may occur at larger distances (> 1000 AU), i.e. 
outside the field-of-view of our observations. 

If the outflow originates from a disk wind, we can obtain 
a constraint on the launch radius from the up per limit on the 
rotatio n velocity of the outflowing gas. Indeed. lAnderson et alJ 
J200l derived a general relation between the poloidal and 
toroidal velocity components of cold magnetocentrifugal wind 
at large distance of the star and the rotation rate of the launching 
surface, independent of the uncertain launching conditions. This 
relation relies on the following facts: 1) The energy and angular 
momentum in the wind are extracted mostly by magnetic fields 
from the rotating disk; 2) The energy extracted is the work done 
by the rotating disk against the magnetic torque responsible for 
the angular momentum extraction; And 3) most of the wind en- 
ergy and angular momentum at observable distances are in the 
measurable kinetic form. In the simple geometry of HH 30, we 
identify the poloidal and toroidal velocities respectively with the 
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radial velocity V ra d and the rotation velocity V rot at a given ra- 
dius r. The rotation rate of the launching surface is the Keplerian 

frequency Oq = ^GM* ^ ~j . lAnderson et alJ ( 120031) relation: 

V ro t = V r 2 ad /(2Qor), can tnen be transformed into 

( VkepM V rot ) ^ 
Haunch = r|2 ^- 

where VkepM = y/GM*/r is the Keplerian velocity at radius r. 
Our upper limit of V rot < 1 kms 1 at r ~200 AU from the flow 
axis indicates ri auncn < 14 AU. 

In this scenario, the optical jet and the outflow would be 
associated with different parts of the disk wind. The jet cor- 
responds to the inner (< 1 AU), densest parts of the wind, 
which has a smalle r opening angle (5° half-opening angle, 
iBacciotti et all 19991) . larger velocities (around 300 km s _1 ), and 
recollimates closer to the disk (100 AU). The CO outflow would 
correspond to the outer part of the wind, ejected at larger radii 
from the star (5-15 AU), and could recollimate at larger dis- 
tances (» 1000 AU). In terms of energetics, these observations 
would suggest that similar momentum fluxes are extracted from 
the disk at these different radii (see Table|5j. 

Several hypothesis could be proposed to explain the appar- 
ent discontinuity between the jet and the outflow, i.e. the lack of 
detected material at intermediate velocities and angles. We first 
note that the outflow cone is not totally empty; in fact, in the 
outflow modelling presented above, adding some material in the 
inner part of the jet allows us to obtain an even better fit to the 
data, to the price however of several new free parameters. One 
possible explanation is a projection effect: The outflow velocity 
field could lie more and more in the plane-of-sky (and thus in- 
distinguishable from the surrounding cloud) when approaching 
the jet axis. The lack of significant detection could also come 
from temperature gradient, that would make the 12 CO J=2-l line 
weaker at smaller opening angles (higher level transitions would 
however emit more strongly, until the gas becomes primarily 
atomic). Alternatively, this lack of observed emission at inter- 
mediate opening angle may be genuine, indicating that the mass 
loading mechanism into the disk wind is significantly more effi- 
cient at radii < 1 AU (jet emission) and between ~5 and 15 AU 
(CO flow emission, see radius estimate above) than in the inter- 
mediate zone ~1 to 5 AU. One could further speculate that the 
region of the disk at the origin of the weak loading level is cooler, 
being e.g. shadowed by an inner puffed-up rim. 

Finally, in this disk wind scenario, the strong difference be- 
tween the northern and southern CO lobes must be related to in- 
trinsic differences in the ejection mechanism due e.g. to physical 
differences between the northern and southern faces of the disk, 
or to a different magnetic configuration. HH30 also presents 
other, second-order, asymmetries. Optical studie s indicate that 
the je t is currently perpendicular to the disk axis dBurrows et al.l 
our measurements however suggest a slight misalignment 
between the flow axis and the jet/disk axis, the inclinations dif- 
fering by 6 + 3°. In the MHD wind model discussed here, this 
suggests the magnetic field does not have a symmetric pattern 
around the disk axis. An additional complexity is the fact that, 
at slightly larger scales (~ 3000 AU), the jet is w obbling around 
its mean direction by ±3° (iBurrows et alj[l996l) . This is an in- 
dication of a precession (or more generally wandering) of the 
ejection direction. However, the data presented in this paper do 
not probe the properties of the molecular outflow at this scale. 



8. Conclusions 

In the HH 30 case, a judicious use of the complementary line 
emission allowed us to derive a number of very robust results: 1) 
The systemic velocity of HH 30 is 7.25 km s ; 2) The disk outer 
radius is 420 AU; 3) The rotation vector of the disk points to- 
ward the North-Eastern jet; 4) The disk is in Keplerian rotation; 
5) The stellar mass is 0.45 M ; 6) A highly asymmetric molec- 
ular outflow originates from the inner parts of the disk; 7) The 
outflow material is mainly located on the thin edges of a cone 
with an opening angle of 30°; 8) The outflow velocity is essen- 
tially radial with a magnitude of 12 km s _1 ; 9) No rotation of the 
outflow is detected. 

The total disk mass is less certain: ~A 10~ 3 M . This depends 
on the assumed dust emissivity. Our measurement indicate a 
rather low value of /?, 0.4, suggestive of significant grain growth. 
In addition, but with still less certainty, the disk appears cold, 
around 12 K at 100 AU, suggesting a significant temperature gra- 
dient between the disk plane and the surface layers. Despite this 
low temperature, there is no measurable depletion of CO com- 
pared to the average abundance in the Taurus cloud, in contrast 
with the general trend for a ll other circumste llar disks (except 
those around HAe stars, see iPietu et al!2006l) . 

The direct determination of the stellar mass improves the 
evolutionary status of HH 30, but this remains hampered by the 
uncertainties on the total luminosity. The most plausible solution 
is that of a ~4 Myr old object, with a total luminosity of order 
0.2 L , but a younger and brighter star cannot be excluded. 

Finally, the origin of the molecular outflow remains unclear, 
with no conclusive arguments to distinguish between an entrain- 
ment mechanism by the optical jet or a disk wind. The one sided 
nature of the outflow remains an unexplained issue in both hy- 
pothesis. 

Acknowledgements. We acknowledge the IRAM staff at Plateau de Bure and 
Grenoble for carrying out the observations. JP thanks S. Cabrit for useful com- 
ments on jet and outflow models and M. J. Welch for providing the 13 CO J=l-0 
spectra cube of the HL Tau environment. We thank both the editor C. Bertout and 
the referee R. Bachiller for theirs comments that helped to improve the paper. 



References 

Anderson, J. M., Li, Z.-Y., Krasnopolsky, R., & Blandford, R. D. 2003, ApJ 

Lett., 590, L107 
Appenzeller, L, Bertout, C, & Stahl, O. 2005, A&A, 434, 1005 
Arce, H. G., Shepherd, D., Gueth, F., et al. 2006, in Protostars and Planets V, 

http://www.ifa.hawaii.edu/UHNAI/ppv.htm 
Bacciotti, R, Eisloffel, J., & Ray, T. P. 1999, A&A, 350, 917 
Bacciotti, R, Ray, T. P., Eisloffel, J., et al. 2003, Ap&SS, 287, 3 
Bachiller, R., Martin-Pintado, J., Tafalla, M., Cernicharo, J., & Lazareff, B. 1990, 

A&A, 231, 174 

Baraffe, I., Chabrier, G., Allard, R, & Hauschildt, P. H. 1998, A&A, 337, 403 
Burrows, C. J., Stapelfeldt, K. R., Watson, A. M., et al. 1996, ApJ, 473, 437 
Cernicharo, J. & Guelin, M. 1987, A&A, 176, 299 

Coffey, D., Bacciotti, R, Woitas, J., Ray, T. P., & Eisloffel, J. 2004, ApJ, 604, 
758 

Cotera, A. S., Whitney, B. A., Young, E., et al. 2001, ApJ, 556, 958 
D'Antona, R & Mazzitelli, I. 1997, Memorie della Societa Astronomica Italiana, 
68, 807 

Dartois, E., Dutrey, A., & Guilloteau, S. 2003, A&A, 399, 773 
Draine, B. T. 2006, ApJ, 636, 1114 

Dutrey, A., Guilloteau, S., & Ho, P. 2006, in Protostars and Planets V, 

http://www.ifa.hawaii.edu/UHNAI/ppv.htm 
Dutrey, A., Guilloteau, S., & Simon, M. 1994, A&A, 286, 149 
Gueth, F. & Guilloteau, S. 1999, A&A, 343, 571 
Gueth, R, Guilloteau, S., & Bachiller, R. 1996, A&A, 307, 891 
Guilloteau, S. & Dutrey, A. 1998, A&A, 339, 467 
Kenyon, S. J., Dobrzycka, D., & Hartmann, L. 1994, AJ, 108, 1872 
Koerner, D. W., Sargent, A. I., & Beckwith, S. V. W. 1993, Icarus, 106, 2 
Konigl, A. & Pudritz, R. E. 2000, Protostars and Planets IV, 759 



16 



Pety et al.: Plateau de Bure Interferometer Observations of the Disk and Outflow of HH 30 



Mannings, V., Koerner, D. W., & Sargent, A. I. 1997, Nat, 388, 555 
Masson, C. R. & Chemin, L. M. 1992, ApJ Lett., 387, L47 
Natta, A., Testi, L., Neri, R., Shepherd, D. S., & Wilner, D. J. 2004, A&A, 416, 
179 

Palla, F. & Stahler, S. W. 1999, ApJ, 525, 772 

Pety, J. 2005, in SF2A-2005: Semaine de l'Astrophysique Francaise, ed. 

F. Casoli, T. Contini, J. M. Hameury, & L. Pagani, 721-+ 
Pietu, V., Dutrey, A., & Guilloteau, S. 2006, in preparation 
Pringle, J. E. 1981, ARAA, 19, 137 

Roddier, C, Roddier, F, Northcott, M. J., Graves, J. E., & Jim, K. 1996, ApJ, 
463, 326 

Shu, F. H., Najita, J. R., Shang, H., & Li, Z.-Y. 2000, Protostars and Planets IV, 
789 

Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593 

Simon, M., Dutrey, A., & Guilloteau, S. 2000, ApJ, 545, 1034 

Stahler, S. W. 1994, in ASP Conf. Ser. 65: Clouds, Cores, and Low Mass Stars, 

ed. D. P. Clemens & R. Barvainis, 330-+ 
Stapelfeldt, K. & The WFPC2 Science Team. 1997, in Science with the VLT 

Interferometer, 395-+ 
Stapelfeldt, K. R. & Padgett, D. L. 1999, in Science with the Atacama Large 

Millimeter Array (ALMA) 
Stapelfeldt, K. R., Watson, A. M., Krist, J. E., et al. 1999, ApJ Lett., 516, L95 
Testi, L., Natta, A., Shepherd, D. S., & Wilner, D. J. 2003, A&A, 403, 323 
Watson, A. M. & Stapelfeldt, K. R. 2004, ApJ, 602, 860 
Welch, W. J., Hartmann, L., Heifer, T., & Briceno, C. 2000, ApJ, 540, 362 
Wood, K., Wolff, M. J., Bjorkman, J. E., & Whitney, B. 2002, ApJ, 564, 887 
Wood, K., Wolk, S. J., Stanek, K. Z., et al. 2000, ApJ Lett., 542, L21 



Pety et al.: Plateau de Bure Interferometer Observations of the Disk and Outflow of HH 30, Online Material p 1 



Online Material 



Pety et al.: Plateau de Bure Interferometer Observations of the Disk and Outflow of HH 30, Online Material p 2 




Fig. 11. Comparison of the channel maps of the '-CO J=2-l observations (left) and our best model (right) of the outflow. The cross 
indicates the disk center as measured in the 1 .3 mm continuum emission. Velocity of the channels are written in km s _1 on the top 
left corner of each channel. Both channel maps share the same color scale and contour spacing (0. 1 Jy/beam). 
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Fig. 12. Comparison of the position-velocity diagrams perpendicular to the jet axis of the 12 CO J=2-l observations (left) and our 
best model (right) of the outflow. The coordinate along the jet axis (5y) are written in arcsecond on the top left corner of each panel. 
All the diagrams share the same color scale and contour spacing (0.1 Jy/beam). 
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Fig. 13. Comparison of the position-velocity diagrams along the jet axis of the 12 CO J=2-l observations (left) and our best model 
(right) of the outflow. The coordinate perpendicular to the jet axis (6x) are written in arcsecond on the top left corner of each panel. 
The horizontal line indicates the position of the outside edge of a cone of 30° half opening-angle. All the diagrams share the same 
color scale and contour spacing (0.1 Jy/beam). 
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Fig. 14. Influence of the different parameters of the outflow model. The left column shows our reference model whose parameters 
are given in Table [5] The each column shows a model for which one and only one parameter (noted on the column top) has been 
varied compared to our reference model. The integrated emission map, the position-velocity diagrams perpendicular and parallel to 
the jet axis and the spectrum integrated over the map are shown from top to bottom. 



Pety et al.: Plateau de Bure Interferometer Observations of the Disk and Outflow of HH 30, Online Material p 5 



V„. = 1 km/a AV = km/a vi a 0.!i" U s = 6" 




S 10 IS 5 10 IS D 5 10 IS OS 10 LS 

VMatLLj (am.B ) VHoeLLy (kra.B ) V#UKLLy (khi.is ) VelasLLy [kjb.b" ) 

Fig. 15. Continuation of Fig. ^] with different varying parameters. 



